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Abstract. The formulation of the tensor virial equations is generalized to unrelaxed con- 
figurations, where virial equilibrium does not coincide with dynamical (or hydrostatic) 
equilibrium. Special classes of homeoidally striated ellipsoids are shown to exhibit similar 
properties as Jacobi ellipsoids, and reduce to Jacobi ellipsoids in the limiting situation of 
homogeneous matter distribution, rigid rotation, and isotropic residual (i.e. other than sys- 
tematic rotation) velocity distribution. Accordingly, the above mentioned density profiles 
are defined as homeoidally striated, Jacobi ellipsoids. Further investigation is devoted to 
the generation of sequences of virial equilibrium configurations. Sequences with constant 
anisotropy parameters are studied with more detail, including both flattened and elon- 
gated, triaxial (in particular, both oblate and prolate, axisymmetric) configurations, and 
the determination of the related bifurcation points. The explicit expression of a number of 
rotation parameters, used in literature, and the ratio of rms rotational to rms residual veloc- 
ity, is also calculated. An application is made to dark matter haloes hosting giant galaxies 
(M w lO 12 iri0), with regard to assigned initial and final configuration, following and 
generalizing to many respects a procedure conceived by Thuan & Gott (1975). The depen- 
dence of the limiting axis ratios, below which no configuration is allowed for the sequence 
under consideration, on the change in mass, total energy, and angular momentum, during 
the evolution, is illustrated in some representative situations. The dependence of the axis 
ratios, £31 and £21, on a parameter, related to the initial conditions of the density perturba- 
tion, is analysed in connection with a few special cases. The same is done for the rotation 
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parameters and the ratio of rms rotational to rms residual velocity. Within the range of the 
rotation parameter, exp 10 (log A =F f"io g A) = exp 10 (— 1.3761 =F 0.500), consistent with 
high-resolution numerical simulations, the shape of dark matter haloes is mainly decided 
by the amount of anisotropy in residual velocity distribution (the above formula uses the 
general notation, exp a (;r) = a x , where the special case, a = e, reduces to the standard 
notation). On the other hand, the contribution of rotation has only a minor effect on the 
meridional plane, and no effect on the equatorial plane, as bifurcation points occur for 
larger values of A. To this respect, dark matter haloes are found to resemble giant elliptical 
galaxies. 

Key words: Cosmology: dark matter - galaxies: evolution - galaxies: formation - galax- 
ies:haloes - galaxies: structure. 

1. Introduction 

Homeoidally striated density profiles make a useful tool in the description of self-gravitating 
fluids on many respects. First, the related potential-energy tensors are expressed as the product 
between a profile factor and a shape factor. Accordingly, a physical parameter which depends 
on the ratio between different components of a same potential-energy tensor, is a function of 
the axis ratios only (Roberts 1962). In addition, the constraints for homeoidally striated ellip- 
soids, in connection with both equilibrium and stability, are left unchanged with respect to the 
special case of homogeneous configurations with same shape (e.g., Binney 1978; Pacheco et 
al. 1989). Second, a criterion for an upper limit to the point of bifurcation from axisymmetric 
to triaxial configurations, is formally identical for collisional and collisionless fluids (Wie- 
gandt 1982a). On the other hand, a condition for dynamical stability depends on the degree of 
anisotropy of peculiar velocity distribution (Wiegandt 1982b). In this framework, additional 
features may be included, such as vorticity and streaming motions (e.g., Pacheco et al. 1986, 
1989; Busarello et al. 1989, 1990; and further references therein). 

For the purpose of the current attempt, our attention shall be restricted to the simplest 
case where both vorticity and streaming motions are not explicitly investigated, and only the 
contribution from systematic rotation is considered separately. Owing to a generalization of 
Dedekind's theorem, to each configuration with no vorticity and assigned angular velocity, 17, 
will correspond an "adjoint" configuration with same boundary and distribution of peculiar 
velocity, but with no rotation and vorticity parameter, Z = ft (Pacheco et al. 1989). Accord- 
ingly, to each transition from an initial to a final, (at least rigidly) rotating configuration, will 
correspond a transition from an initial to a final, vorticating, adjoint configuration. 

The choice of homeoidally striated density profiles as a valid approximation to the descrip- 
tion of self-gravitating fluids, is motivated by two main reasons. In primis, most astronomical 

©WILEY- VCH Verlag Berlin GmbH, 13086 Berlin, Germany 2004 0044-6337/04/32201-0223 $ 17.50+.50/O 



R. Caimmi and C. Marmo: homeoidally striated density profiles 



225 



bodies exhibit ellipsoidal-like shapes and isopycnic surfaces. Secundo, homeoidally striated 
ellipsoids approximate to a first extent collisional and collisionless, self-gravitating fluids (e.g., 
Vandervoort 1980; Vandervoort & Welty 1981; Lai et al. 1993), at least with regard to global 
properties, involving averages on the whole mass distribution. On the other hand, some caution 
must be used in dealing with local properties, involving values in situ. For instance, velocity 
dispersion values measured in the central region of a galaxy, can be significantly different for 
structurally identical systems (and so characterized by identical, global velocity dispersions), 
due to different orbital structures. 

Though classical investigations have been restricted to figures of dynamical (or hydro- 
static) equilibrium such as MacLaurin spheroids and Jacobi ellipsoids (e.g., Jeans 1929, 
Chaps. VIII-IX; Chandrasekhar 1966, Chaps. 5-8), still the current attempt may be general- 
ized to figures of virial equilibrium. Let us define virial equilibrium as characterized by the 
validity of the virial theorem, and relaxed and unrelaxed configurations as systems where 
virial equilibrium does and does not coincide, respectively, with dynamical (or hydrostatic) 
equilibrium. 

In the special case of relaxed configurations, the virial theorem is usually formulated by 
setting the sum of the potential energy, and twice the kinetic energy, equal to zero. On the 
other hand, the above mentioned relation no longer holds for unrelaxed configurations, and 
a more general formulation of the virial equations is needed. To this aim, our attention shall 
be limited to homeoidally striated density profiles, where the rotational velocity field is sub- 
jected to a number of necessary restrictions. Then the related sequences of virial equilibrium 
configurations with constant rotation and anisotropy parameters, can be studied in detail. 

An extension of the formulation of the virial equations to unrelaxed configurations, and the 
investigation of sequences where both rotation and anisotropy parameters remain unchanged, 
make the aim of the current paper. As an application of the theory, a procedure to establish 
mutual connection between physical parameters related to assigned initial and final state of a 
density perturbation, may be generalized to many respects. 

About a quarter of century ago, Thuan & Gott (1975) idealized elliptical galaxies as 
MacLaurin spheroids, resulting from virialization after cosmological expansion and subse- 
quent collapse and relaxation of their parent density perturbations. Unexpectively, the above 
mentioned model provided an answer to the question, why there are no elliptical galaxies more 
flattened than E7. The method is based on the knowledge of an initial configuration, assumed 
to be at turnaround, and a final configuration, assumed to be a MacLaurin spheroid. A gener- 
alization to a somewhat more arbitrary initial and final configuration, and to other respects, is 
expected to be useful for e.g., comparison with (i) results of numerical simulations, and (ii) 
observations. 

More specifically, the problem may be formulated as follows. Let the physical parameters 
of a density perturbation be specified at a selected, initial and final configuration, at redshift 
Zi and Zf, respectively, provided z rec > Zi > Zf > z vir , where the indices, rec and vir, 
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denote recombination and virialization, respectively. Let the density profile be assigned, and 
the isopycnic surfaces, i.e. the surfaces of equal density, be similar and similarly placed ellip- 
soids. Let the rotational velocity field, and the amount of anisotropy in residual (in particular, 
peculiar) velocity distribution, also be specified. Under the above mentioned assumptions, de- 
termine the shape of the final configuration, and the related physical parameters. 

The generalization of Thuan & Gott (1975) procedure to a special class of homeoidally 
striated ellipsoids, together with an application to dark matter haloes hosting giant galaxies 
(M « lO 12 ni0), makes a specific application of the theory outlined above. It shall be analysed 
in the present attempt. 

The current paper is organized as follows. The general theory of homeoidally striated den- 
sity profiles, and the results of interest to the aim of this investigation, are outlined in Sect.|2] 
The properties of a special class of configurations, defined as homeoidally striated Jacobi el- 
lipsoids, are studied in Sect.[3] where, in particular, the explicit expressions of some rotation 
parameters, known in literature, are derived. Then sequences of virial equilibrium configu- 
rations are determined, where both rotation and anisotropy parameters remain unchanged. A 
generalization of the procedure followed by Thuan & Gott (1975) to homeoidally striated den- 
sity profiles is performed in Sect.|4l and a number of special situations are also investigated. 
Furthermore, an application is made to dark matter haloes hosting giant galaxies, and the re- 
sults are presented and discussed. Some concluding remarks are drawn in Sect. [5] and a few 
arguments are treated with more detail in the Appendix. 

2. General theory 

A general theory for homeoidally striated density profiles has been developed in earlier ap- 
proaches (Roberts 1962; Caimmi 1993a; Caimmi & Marmo 2003, hereafter quoted as CM03), 
and an interested reader is addressed therein for deeper insight. What is relevant for the current 
investigation, shall be mentioned and further developed here. 
The isopycnic surfaces are defined by the following law: 



where po, ro are a scaling radius and a scaling density, respectively, related to a reference 
isopycnic surface, and S, R, correspond to the boundary. By "radius" we intend here the radial 
coordinate of a generic point on a generic isopycnic surface. The radius changes - in general - 
passing from one point to one other, along a fixed isopycnic surface, and it may coincide - in 
particular - with a semiaxis of the related ellipsoid. 

Let Si and £2 be generic isopycnic surfaces belonging to a homeoidally striated density 
profile, and n, ri, generic radii along the same direction. It can be shown (e.g., CM03) that 
the ratio, r\ jr-x, is independent of direction. The above result holds, in particular, in the special 



P = pof(0 ; /(i) = i ; po=pQ) ; 
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case where an isopycnic surface coincides with the reference isopycnic surface, according to 
Eq. dlb> . In other words the scaled radius, £, is defined as the ratio between the radius of 
the isopycnic surface under consideration, along a generic direction, and the radius of the 
reference isopycnic surface, along the same direction. 

The mass, the inertia tensor, and the self-energy tensor, are: 

M = v mas M Q ■ (2) 

Ipq = SpqVinrMap \ (3) 

GM 2 

(Esel)pq — v se i(B se i) pq ; (4) 

ax 

3 

(Bsei)pq = 5 pq tp2e p zA p ; B se i —^^t S 2e si A s ; (5) 

s=l 

where 8 pq is the Kronecker symbol; G is the constant of gravitation; 
profile factors i.e. depend only on the density profile via the scaled radius, S; aj., 03, are 
semiaxes; e pq = a p /a q are axis ratios; A\, A%, A3, are shape factors i.e. depend only on the 
axis ratios; and Mo is the mass of a homogeneous ellipsoid with same density and boundary 
as the reference isopycnic surface: 

Mo = — /5 aoiao2ao3 I (6) 
where a i, a 02 , a 03 , are the semiaxes of the ellipsoid bounded by the reference isopycnic 
surface. Owing to a different normalization of the physical quantities, the profile parameters, v, 
have a different expression with respect to CM03. For further details, see Appendix A. Values 
of shape factors related to limiting configurations (round, oblate, prolate, flat, and oblong) are 
derived in Appendix B. 

In dealing with angular momentum and rotational energy, the preservation of (triaxial) el- 
lipsoidal shape imposes severe constraints on the rotational velocity field. Leaving an exhaus- 
tive investigation to more refined approaches, our attention shall be limited here to a restricted 
number of special cases, namely: (i) rigid rotation about a principal axis, and (ii) differential 
rotation about a symmetry axis. By "rotation" it is intended, of course, circular rotation. 

In the special case of rigid rotation about a principal axis, let it be x 3 , the angular- 
momentum vector and the rotational-energy tensor are: 

J s = S s3 (l pp + I qq )Q s = S s3 u mr Ma p (l + e qp )(v rot ) p ; p ^ q ^ s ; (7) 

(Erot)pq = -^Ipq&sZ&l = ^^Pl^ ~~ &ps)Vinr M a p Q 2 
1 2 

= 2<^3<Wl ~ S ps)v lnr M [(v rot ) p ] ; Pt^s] q ^ s ; (8) 
where is the angular velocity related to the axis, x 3 , and (v ro t) p is the rotational velocity at 
the end of the semiaxis, a p , p = 1,2. 

The module of the above vector and the trace of the above tensor, which make the related 
angular momentum and rotational energy, respectively, read: 

J = u lnr Ma p (l + e qp )fl = v inr Ma p (l + e qp )(v rot ) p ; (9) 
Erot = ~v inr Ma 2 p {l + e 2 qp )n 2 = ~v inr M{l + e 2 qp ) [{v rot ) p } 2 ; (10) 
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for further details, see CM03. 

In the special case of differential rotation about a symmetry axis, let it be X3, our attention 
shall be restricted to rotational (with respect to x 3 axis) velocity distributions which obey the 
law: 

Vrot{r,0) = Vrot(a' p ,0) 

v rot (R,9) v ro t{a p ,0) 

or equivalently, to angular (with respect to x 3 axis) velocity distributions which obey the law: 

"(r,g) = "K,0) . (m 

n{R,e) n(op,o) ' 

where (r, 6), (R, 9), represent a point on a generic isopycnic surface and on the boundary, 
respectively, along a fixed radial direction, and (a' , 0), (a p , 0), p = 1, 2, represent the end of 
the corresponding equatorial semiaxis. It is worth noticing that, in particular, either rotational 
or angular velocity are allowed to be constant everywhere. 

The equivalence of Eqs. (lilt and fl!2i is owing to the similarity between isopycnic sur- 
faces, which translates into the equivalent relations: 

^ = a -; ^= a -; (13) 

R a ro ao 

as outlined above in connection with ellipsoidal configurations. 

The angular-momentum vector and the rotational-energy tensor are: 

J s = 2S s3 7] anm v anm Mav rot ; (14) 

{E ro t) P q = S pq (l - Sp^rjrotVrotMVrot | (15) 

where v rot is the rotational velocity at the end of the equatorial semiaxis, a; rj anm , r) ro t, depend 
on the boundary and on the distribution of angular velocity therein, and then may be conceived 
as shape factors; on the other hand, v anm , v rot , depend on the density profile and on the radial 
distribution of angular velocity, and then are genuine profile factors. 

The module of the above vector and the trace of the above tensor, which make the related 
angular momentum and rotational energy, respectively, read: 

J = 2r) anm v anm Mav rot ; (16) 

E rot = 1r\ TO tv TO tMv1 ot ; (17) 

for further details, see CM03. 

To save space, let us merge Eqs. and d!4l >: (|SJl and dl5l i: as: 

Js = S s3 i] anm iy anm Ma p (l + e 2 qp ){vrot) P \ p ^ q ^ s ; (18) 

(Erot)pq = S pq (l - 8 p z)rj rot v ro tM [(tvot)p] 2 ; (19) 
and Eqs. © and QTjJ; iQOj and (Q3; as: 

J = r} anm v anm M (1 + £2i) a i(*Vot)i ; (20) 

Erot = r] rot v rot M{\ + e% x ) [(v rot )i} 2 ; (21) 

respectively, where rj anm i/ anm = Vi nr and r\ ro tv ro t — v. mr /2 in the special case of rigid 
rotation; for further details, see CM03. 
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An expression of the rotational energy, E rot , as a function of the angular momentum, J, 
is obtained by the combination of Eqs. (I20i and OIL as: 

J_ 2 

Ma{ 

V T ot 



Erot — n T 2^ rarn ^ ram ' (22a) 



(22b) 



i e 2 

\^Tarn)pq ^pq\^- ^P^~^2 7l i 2 \2 5 (22c) 



?7rot 1 

Vanm 1 + e 21 

and the ratio of rotational to self-potential energy, taken with the positive sign, £ rot , is obtained 
by the combination of Eqs. @, extended to the trace of the related tensor, and \22\ . as: 

Erot V r am E r am , * 

h ; (23a) 



E se l Vsel B 

J 2 



(23b) 



GAP ai ' 

where both f rot (e.g., Ostriker & Peebles 1973) and h (e.g., Caimmi 1980) may be conceived 
as rotation parameters. 

At this stage, we have all the ingredients for generalizing the usual formulation of the 
tensor virial equations (e.g., Chandrasekhar 1961, Chap. 13, §117; Brosche 1970; Brosche 
et al. 1983), to include unrelaxed configurations, in particular homeoidally striated density 
profiles. 

3. The tensor virial equations for homeoidally striated density profiles 
3.1. Basic theory 

As outlined in Sect. 1, homeoidally striated density profiles make a useful tool in the descrip- 
tion of self-gravitating fluids on many respects. Though classical investigations have been 
restricted to figures of dynamical (or hydrostatic) equilibrium such as MacLaurin spheroids 
and Jacobi ellipsoids (e.g., Jeans 1929, Chaps. VIII-IX; Chandrasekhar 1966, Chaps. 5-8), still 
the current attempt may be generalized to figures of virial equilibrium. To this respect, it is 
worth remembering that the virial theorem holds for values of parameters averaged over a suf- 
ficiently long time, t T, where T is a characteristic period of the system (e.g., Landau & 
Lifchitz 1966, Chap. II, § 10). 

Let us define virial equilibrium as characterized by the validity of the virial theorem, and 
relaxed and unrelaxed configurations as systems where virial equilibrium does and does not 
coincide, respectively, with dynamical (or hydrostatic) equilibrium. For instance, a compress- 
ible homogeneous sphere undergoing coherent oscillations, is in virial equilibrium: then the 
generic configuration is unrelaxed, even if a special one exists, equal in shape to its relaxed 
counterpart. 

In the special case of relaxed configurations, the virial theorem is usually formulated by 
setting the sum of the potential energy, and twice the kinetic energy, equal to zero. On the 
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other hand, the above mentioned relation no longer holds for unrelaxed configurations, and a 
more general formulation of the virial equations is needed. To this aim, our attention shall be 
limited to homeoidally striated density profiles, where the rotational velocity field is subjected 
to the restrictions mentioned in Sect. [2] 

With regard to unrelaxed configurations, let us define the residual kinetic-energy tensor, 
(E res ) pq , as owing to any motion other than systematic rotation e.g., random motions, stream- 
ing motions, radial motions; and the effective, residual kinetic-energy tensor, (E res ) pq , as the 
right amount needed for the configuration of interest to be relaxed i.e. the validity of the tensor 
virial equations: 

(Esel)pq + 2(E ro t)pq + 2(E res ) pq — ; (24) 

then the components of the effective anisotropy tensor, defined as: 

C Pq = {Eres)pq ; p= 1,2,3; 9 = 1,2,3; (25a) 

3 

Cpp = 1 ; < C PP < 1 ; C Pq = ; p # q ; (25b) 

P =i 

may be conceived as effective anisotropy parameters (with respect to the effective residual 
kinetic-energy tensor). Similarly, the components of the generalized anisotropy tensor, defined 
as: 



Cpq = {Eres)pq ■ p= 1,2,3; g = l,2,3; (26a) 

±L/ res 

3 w 



E 

Cpp = -rr^- = C ; < C PP < C ; ( pq = ; p ^ q ; (26b) 

p=l 



may be conceived as generalized anisotropy parameters (with respect to the generalized resid- 
ual kinetic-energy tensor). 

The scalar, £, may be thought of as a virial index, where £ = 1 corresponds to null virial, 
(E se i) pq + 2(E rot ) P q + 2(E res ) pq = (which doe not necessarily imply a relaxed configura- 
tion), ( > 1 to negative virial, and £ < 1 to positive virial. 

The combination of Eqs. d25l > and i26\ yields: 

Cpp = ^ ; P = 1,2,3 ; (27) 

and the specification of the effective and generalized anisotropy parameters, ( pp and ( pp , and 
the residual kinetic energy, E res , implies the specification of the effective, residual kinetic- 
energy tensor, (E res ) pq . 

In the special case of relaxed configurations, the effective, residual and actual kinetic- 
energy tensors coincide, (E res ) pq = (E res ) pq , i.e. £ = 1. Then the limiting situation of 
relaxed configurations may directly be obtained from the results of the current Section, in the 
limit C ->■ 1. 

Using Eqs. @, (0, (II 91 . (I22> . and (I26> . Eq. 1241 takes the more explicit expression: 

GM 2 J 2 

b'sel (E sel)pp ~t~ 2v TaTa — — n (-^ramjpp 2C >p pE res — , (28) 

a\ ivia^ 



R. Caimmi and C. Marmo: homeoidally striated density profiles 



231 



leaving aside the trivial case of non-diagonal terms, where they reduce to identities, = 0. 

The special choice, p — 3, reads: 

GM 2 

- Vsel {Bsel ) 3 3 + K^E res = ; (29) 

which allows the following expression of the residual kinetic energy: 

E res = -^-r^~ (B se i) 33 ; (30) 

2 C33 ai 

as first outlined in dealing with MacLaurin spheroids (Brosche 1970). 



The combination of Eqs. J28i and d30i yields 



Vse.1 



GM 



2 r 



ax 



(B se l)qq ~ y^-{B se l) 33 
(,33 



J 2 



;(B ram ) qq = Q ; 9=1,2 ; (31) 



Ma{ 

leaving aside the trivial case, p = 3, where it reduces to an identity, = 0. A necessary 
condition for satisfiyng Eq. J3 II . or Q4> . using Eqs. 0, may be written as: 

(qq_ , (Bsel)qq _ A q _ _ 
(,33 \&sel )33 e 3q A 3 

which is the natural extension of its counterpart related to axisymmetric, relaxed configura- 
tions (Wiegandtl982a-b). 

3.2. Rotation parameters 

The amount of rotation parameter, h, needed to satisfy Eqs. d28l >. can be determined via 

Eqs. J22bt and J23bi by use of Eq. ( I3U . The result is: 

1 v se i (B se i) qq 



2 Vram {,B rarn ) qq 



^ (qq {Bsel)33 



C33 {B se l) 



qq 



9 = 1,2; (33) 

and the related value of the rotation parameter, £ ro t, owing to Eqs. (j53aj and j33L takes the 
explicit expression: 

c 1 -B ram {B se l) qq [_ | (-Bsei )33 



2 -Bsei (-B r 



C33 {Bseljq 



9 = 1,2; (34) 
where £ gg = (33 = 1/3 for isotropic, residual velocity distributions (i.e. related to the residual 
kinetic energy); e 9 2 = £gi = 1 for spheroidal configurations; (l+e|i)-B ram = X/2,v ram = 5, 
for rigidly rotating configurations; i/ se ; = 3/10 for homogeneous configurations. 
The validity of Eq. (I33> . or equivalently Eq. (I34L implies the following relation: 

C33 [(-Bsei) 11 (A-am) 22 — {B S el)22{B ra m)ll\ = (^sei)33 [Cll (-E?r am) 22 ~ C22 {B r am)ll] 1 (35) 

which links the axis ratios, £21 and £31, with the anisotropy parameters, (n and £22, owing 
to Eq. i26\ . For fixed values of the (generalized or effective) anisotropy parameters, Eq. i35\ 
describes the change of an axis ratio as a function of the other one, along the related sequence 
of virial equilibrium configurations. 

The rotation parameters, h and £ ro t, expressed by Eqs. J33i and 1341 . respectively, exhibit 
a substantial difference: while the former depends on profile factors, v se i and v ram , the lat- 
ter does not. Then the value of £ rot remains unchanged for homeoidally striated ellipsoids 
with fixed axis ratios and anisotropy parameters. With regard to h, the special case of Jacobi 
ellipsoids reads (e.g., Chandrasekhar 1962, but the notation is different therein): 

h = ^ 1+ 2 2l) l(B S el)qq-(B S elh3} i 9 = 1,2; (36) 

t>u e ql 
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and the special case of MacLaurin spheroids corresponds to £21 = 1, (B ram ) qq = 1/8, owing 
to Eq. (I22c> and rigid rotation. For further details on MacLaurin spheroids and Jacobi ellip- 
soids see e.g., Jeans (1929, Chap. VIII); Chandrasekhar (1969, Chaps. 5-6); Caimmi (1996a). 

At this stage, let us define the rotation parameter: 

ft 2 

" = ; (37) 

where ft = f2(ai, 0) is the angular velocity at the end of the major equatorial semiaxis, de- 
noted as m, and p is the mean density of the ellipsoid: 

47r ai<22a3 

the above definition of the rotation parameter, v, makes a generalization of some special cases 
mentioned in the literature (e.g., Jeans 1929, Chap. IX, §232; Chandrasekhar & Leboviz 1962). 
The combination of Eqs. J20l >. J23bi . d37i . and J38I yields: 

h = -n 2 v 2 ( 1 + £ 2l) 2 . (m 

^ £21^31 

which links the rotation parameters, h and v. 

The combination of Eqs. ( I22bl i. (22c\ . d33t . and d39l allows the explicit expression: 

Cqq {B se l)33 



1 v sei 

«' Vrot^rot 



1 



(7=1,2; (40) 



C33 (B se l)qq_ 

which, in the special case of rigidly rotating, homogeneous configurations with isotropic pe- 
culiar velocity distribution, reduces to a known relation for Jacobi ellipsoids and, with the 
additional demand of axial symmetry, to a known relation for MacLaurin spheroids (e.g., 
Jeans 1929, Chap. VIII, §§189-193; Chandrasekhar 1969, Chap. 5, §32, Chap. 6, §39; Caimmi 
1996a). 

In conclusion, Jacobi ellipsoids and MacLaurin spheroids may be conceived as limiting 
cases within the class of homeoidally striated ellipsoids under discussion i.e. characterized 
by rotational or angular velocity distribution expressed by Eqs. Jl Q or dl2i . respectively. For 
this reason, we define homeoidally striated ellipsoids or spheroids of the kind considered, as 
homeoidally striated Jacobi ellipsoids or MacLaurin spheroids, respectively. 

On the other hand, one must keep in mind that Jacobi ellipsoids and MacLaurin spheroids 
are relaxed configurations i.e. in dynamical equilibrium (e.g., Jeans 1929, Chap. VIII, § § 1 89- 
193; Chandrasekhar 1969, Chap. 5, §32, Chap. 6, §39), while their homeoidally striated coun- 
terparts are unrelaxed configurations i.e. in virial equilibrium but, at least in general, not in 
dynamical (or hydrostatic) equilibrium (e.g., Tassoul 1978, Chap. 4, Sect. 4.3; Chambat 1994; 
Caimmi 1996a,b). In addition, homeoidally striated ellipsoids approximate to a first extent 
self-gravitating fluids, at least in connection with typical values of physical parameters, aver- 
aged over the whole volume (e.g., Vandervoort 1980; Vandervoort & Welty 1981; Lai et al. 
1993). 
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3.3. Isotropic and anisotropic configurations 



The ratios between anisotropy parameters can be obtained via Eq. ( I33l l. The result is: 
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and the combination of Eqs. (1261 and d41l yields: 
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which provides an alternative expression to Eqs. (I33L ( I34l >. and (1401 . as: 
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(41) 
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where the rotation parameter, h, is expressed as a function of profile parameters, shape param- 
eters, and the effective anisotropy parameter, £33 = C33/C I n addition, Eqs. J44l >. ( 145 > . and 
( I46> . have the same formal expression as in the special case of axisymmetric configurations 
(e.g., Caimmi 1996b). 

The combination of Eqs. WW and d43t yields: 
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in the limiting case of a vanishing rotation, h — > 0, Eqs. J47i reduce to: 
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p= 1,2,3; (47) 



(48) 



where flattening and/or elongation are owing to anisotropic, residual velocity distribution (i.e. 
related to the residual kinetic energy) only. 

In the special case of isotropic configurations, Q pp = 1/3, p = 1, 2, 3, Eq. J44i reduces to: 
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(49) 



and an equivalent expression may be obtained from Eq. ( I33l >. The combination of Eqs. (I44t . 
and j49t yields: 



h = hi 



1 V se l 3C33 — 1 (B S el)s3 
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33 
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(50) 



which shows the effect of polar (with respect to equatorial) anisotropy for a fixed shape. 

Values (33 < 1/3 make h < hi SO i.e. less rotation is needed for a fixed flattening, as a 
contribution from anisotropy arises. The anisotropy parameter cannot be lower than a critical 
value, which corresponds to h = 0; if otherwise, the shape under consideration would not 
be related to a configuration of virial equilibrium. On the contrary, values £ 33 > 1/3 make 
h > hiso i-e. more rotation is needed for a fixed flattening, as a contribution from anisotropy 
arises, but in the opposite sense with respect to the former case. 
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Following a similar procedure, Eq. J46> reads: 
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where the effect of anisotropy occurs as discussed above. 

Owing to Eqs. @, d24l. and d26l . the total energy of a homeoidally striated, Jacobi ellip 
soid is: 
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and the combination of Eqs. J23bi . ( I33> . and (1521 yields the explicit expression of the rotation 
parameter, A (e.g., Peebles 1969): 
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an equivalent expression may be obtained using Eq. (1441 instead of J33i : the result is: 
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which is a generalization of an expression, previously calculated for homeoidally striated Ja 
cobi ellipsoids, where £11 = £22 (Caimmi 1993b). 

Let us define, using Eq. j30L a typical rms residual velocity, < v 2 es >= 2E res /M, as: 
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and define, using Eqs. (|SaJ and ( 13 II . a typical square rotational velocity, < v 

2E rot /M, as: 
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then the combination of Eqs. d55t and i5bt yields: 
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which for anisotropic, spheroidal configurations rotating around their symmetry axis, attains 
a different expression with respect to e.g., Binney & Tremaine (1987, Chap. 4, § 3), owing to 
a different definition of < v% es > therein. The rms velocity ratio, Xv, may also be considered 
as a rotation parameter. 

However, some caution has to be used in the interpretation of these results. According 
to Eq. d57i the rms velocity ratio, Xv< depends only on the axis ratios, £21 and £31, and the 
generalized anisotropy parameter, (33 (e.g., Binney 1976; Binney & Tremaine 1987, Chap. 4, 
§ 3). Then homeoidally striated ellipsoids with different density profiles and same generalized 
anisotropy parameter and axis ratios, have the same rms velocity ratio. On the other hand, it is 
not true for self-gravitating fluids, at least in the limit of isotropic, peculiar velocity distribu- 
tions, where the rms velocity ratio tends to zero as the density profile tends to the Roche limit, 
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i.e. a mass point surrounded by a massless atmosphere (Caimmi 1980; 1983). This discrep- 
ancy is probably owing to the fact, that (for an assigned generalized anisotropy parameter) the 
isopycnic surfaces are (taken to be) rotationally distorted at the same extent for homeoidally 
striated ellipsoids, while rotational distorsion increases passing from the centre to the bound- 
ary for self-gravitating fluids. 

3.4. Sequences of homeoidally striated Jacobi ellipsoids 

Sequences of homeoidally striated Jacobi ellipsoids may be generated using Eqs. (I25l > and J35L 
which make a system of three equations involving five unknowns: the effective anisotropy 
parameters, ( pp , and the axis ratios, £21 and £31. If effective anisotropy parameters are kept 
constant along the sequence, then Eq. d35l reduces to a relation involving axis ratios only, and 
the rotation parameter, h, may be determined using Eq. (1401 . If, on the other hand, the rotation 
parameter is kept constant along the sequence, then Eqs. d25i and (I33> make a system of three 
equations in five unknowns, and the allowed values of the effective anisotropy parameters, or 
other rotation parameters, are represented by a surface, z = z(e-2i, £31). 

In the special case of Jacobi ellipsoids, incompressibility acts as an isotropic, residual 
velocity distribution and ellipticities are owing to rotation only, and the related sequence is 
characterized by ( pp = 1/3, p = 1, 2, 3 (e.g., Jeans 1929, Chap. VIII; Chandrasekhar 1966, 
Chaps. 5-8). It is true also for gaseous, relaxed, homeoidally striated Jacobi ellipsoids, where 
the collisional nature of the fluid (in absence of systematic motions different from rotation) 
necessarily yields an isotropic residual velocity distribution, which coincides with Maxwell's 
distribution (e.g., Jeans 1929, Chap. IX; James 1964). On the contrary, an anisotropic resid- 
ual velocity distribution may occur in presence of systematic motions different from rotation 
within collisionless fluids (e.g., Wiegandt 1982 a,b), such as stars or non baryonic dark matter. 

Generally speaking, there is no apparent criterion for the selection of a specified sequence 
of homeoidally striated Jacobi ellipsoids, instead of one other. To this respect, only the nature 
of the special problem under consideration could provide further information. For this reason, 
our attention shall be limited to sequences with constant effective anisotropy parameters, £ pp , 
starting from a nonrotating, axisymmetric configuration. Accordingly, (n = C22, via Eq. J48b . 
The ending point is related to an oblong configuration, unless centrifugal support at the ends 
of the major equatorial axis occurs earlier. 

The equatorial axis ratio, £21, as a function of the polar axis ratio, £31, for sequences 
of density profiles of the kind considered, is plotted in Fig.^ with regard to initial, non- 
rotating (h = 0), spheroidal boundaries characterized by polar axis ratios within the range 
0-1 < (e3i)h=o < 10. The bifurcation from axisymmetric to triaxial configurations is marked 
by a filled circle on each curve. The bifurcation from a round boundary takes place in a se- 
quence where (£31)^=0 = 1.97698. Each sequence is characterized by the initial, nonrotating 
boundary, i.e. the effective anisotropy parameters, £ pp , but it is independent of the special 
choice of density profile and (systematic rotation) velocity profile. 
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E 31 



Fig. 1. The equatorial axis ratio, €21, as a function of the polar axis ratio, 631, for sequences 
of density profiles in virial equilibrium starting from nonrotating (h = 0), spheroidal bound- 
aries, with polar axis ratios (631)^=0 = 0.1 i (oblate) and (eai)h=o = 1/(0.1 i) (prolate), 
1 < i < 10. The bifurcation from axisymmetric to triaxial configurations is marked by a 
filled circle on each curve. The bifurcation from a round boundary takes place in a sequence 
where (631)^=0 = 1.97698. A pentagon skeleton on each curve marks the upper limit of polar 
axis ratio, for which a necessary condition for the occurrence of centrifugal support at the 
ends of the equatorial major axis, is first satisfied with regard to a NFW density profile (see 
Subsect. 4.5) in rigid rotation. 



The same holds for the rotation parameter, £ ro t, defined by Eq. ( 145 \ . With regard to the 
remaining rotation parameters, let us define their normalized counterparts, via Eqs. (1441 . i46i . 
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Fig. 2. The rotation parameters, £ ro u h-N, -V/v, and vn (from top left in clockwise sense), as 
a function of the polar axis ratio, 631, for sequences of density profiles in virial equilibrium, 
starting from nonrotating (h = 0), spheroidal boundaries, with polar axis ratio (e3i)/i=q = 
0.1 i (oblate) and (e3i)h=o = 1/(0-1 *) (prolate), 1 < % < 10. Other captions as in Fig.^ 
All the curves end at the oblong configuration, where the axis ratios tend to zero. The related 
values of the rotation parameters, £ rot and vn, tend to 0.5 and 0, respectively, while and 
A jv tend to infinity. 

2 Xv 

x N = — ; (6i) 

which are also independent of both density profile and (systematic rotation) velocity profile. 

The rotation parameters, £ ro t, /ijv, ^n, and vn (from top left in clockwise sense), as a 
function of the axis ratios, €31 and €21, are plotted in Figs.|2]and[3] respectively, for sequences 
of density profiles in virial equilibrium starting from nonrotating (h = 0), spheroidal bound- 
aries, with polar axis ratios (£31)^=0 = 0.1 i (oblate) and (£31)^=0 = 1/(0.1 i) (prolate), 
1 < i < 10. Each curve is characterized by the effective anisotropy parameters, ( pp , but is 
independent of both density profile and (systematic rotation) velocity profile. 

The rotation parameter, \n, as a function of the axis ratios, £31 and £21, is plotted in 
Fig. for sequences of density profiles in virial equilibrium starting from nonrotating (h = 
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Fig. 3. The rotation parameters, £ ro u ^n, ^n, and vn (from top left in clockwise sense), as a 
function of the equatorial axis ratio, £21, for sequences of density profiles in virial equilibrium, 
starting from nonrotating (h = 0), spheroidal boundaries, with polar axis ratio (£31)^=0 = 
0.1 i (oblate) and (esi)h=o = 1/(0-1 (prolate), 1 < % < 10. Other captions as in Fig.^ 
All the curves end at the oblong configuration, where the axis ratios tend to zero. The related 
values of the rotation parameters, £ rot and Vjy, tend to 0.5 and 0, respectively, while and 
Aat tend to infinity. The coordinates of the nonrotating configurations are (1,0) in each panel. 

0), spheroidal boundaries, with polar axis ratios (£31)^=0 = 0.1 i (oblate) and (£31)^=0 = 
1/(0.1 i) (prolate), 1 < i < 10. 

Axis ratios and normalized rotation parameters related to virial equilibrium configurations 
belonging to three typical sequences, are listed in Tabs.^0 and ??, where the starting nonro- 
tating configuration is oblate (anisotropy induces flattening), spherical (ispotropic case), and 
prolate (anisotropy induces elongation), respectively. 

3.5. A necessary condition for centrifugal support 

With regard to Jacobi ellipsoids, rotational equilibrium i.e. balance between gravitational and 
centrifugal force, first occurs at the end of the major equatorial semiaxis, for fiat or oblong 
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Fig. 4. The rotation parameter, \n, as a function of the polar axis ratio, £31 (upper panel) 
and the equatorial axis ratio, £21 (lower panel), for sequences of density profiles in virial 
equilibrium starting from nonrotating (h = 0), spheroidal boundaries, with polar axis ratios 
(£31)71=0 = 0.1 i (oblate) and (£3i)/ l= o = 1/(0.1 i) (prolate), 1 < i < 10. Other captions 
as in Fig.^ All the curves end at the oblong configuration, where the axis ratios tend to zero 
and the rotation parameter, xn, tends to infinity. The coordinates of the starting configuration 
related to all the sequences are (1, 0) in the lower panel. 

configurations (e.g., Caimmi 1980). On the other hand, homeoidally striated Jacobi ellipsoids 
show rotational equilibrium for less flattened or less elongated configurations. The calculation 
of the gravitational force induced by density profiles of the kind considered, involves numeri- 
cal integrations (e.g., Chandrasekhar 1969, Chap. 3, §20) and is outside the aim of the current 
attempt. 

Our attention shall be restricted here to two limiting situations, for which the calculation 
of the gravitational force induced by the mass distribution is simple, and the corresponding 
range includes the gravitational force induced by the system under investigation. In addition, 
our considerations will be restricted to the boundary, which is the case of physical interest. On 
the other hand, the procedure may be extended to a generic, isopycnic surface. 
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0.0385384 
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0.296287 
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0.0382934 
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0.617888 


0.286432 
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0.275970 


0.75 


0.140272 


0.661862 


0.0374076 


0.999672 


0.624450 


0.264850 
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0.141925 


0.659686 


0.0367259 


1.01438 


0.629567 


0.253015 


0.65 


0.144066 


0.661701 


0.0358531 


1.03357 


0.636204 


0.240404 


0.60 


0.146785 


0.668618 


0.0347618 
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0.644647 


0.226950 


0.55 
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0.0334219 
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0.701221 
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0.668457 
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0.45 


0.159634 


0.729873 


0.0298619 


1.17795 


0.684842 


0.180748 


0.40 


0.166054 


0.769728 


0.0275717 


1.24013 


0.705158 


0.163111 


0.35 


0.173953 


0.824243 


0.0248968 


1.31896 


0.730426 


0.144198 


0.30 


0.183700 


0.898594 


0.0218121 


1.41997 


0.762088 


0.123907 


0.25 


0.195803 


1.00092 


0.0183099 


1.55171 


0.802290 


0.102147 


0.20 


0.211000 


1.14492 


0.0144169 


1.72824 


0.854462 


0.0788405 


0.15 


0.230466 


1.35619 


0.0102269 


1.97547 


0.924692 


0.0539566 


0.10 


0.256360 


1.69132 


0.00596398 


2.35031 


1.02577 


0.0275656 


0.05 


0.294098 


2.32663 


0.00212718 


3.02987 


1.19513 


0.0000000 


0.00 


0.500000 


+oo 


0.00000000 


+oo 


+oo 



Table 1. Values of axis ratios, £31 and £ 2 i, and normalized rotation parameters, £ rot , Jin, Vn, 
-V/v, Xn, related to a sequence of virial equilibrium configurations with anisotropy parameters, 
Cn = C22 = 0.391002 and C33 = 0.217996. The polar axis ratio of the initial nonrotating 
configuration is e 3 i = 1/2. Values at the bifurcation point are listed on the boldface line. 

Given a homeoidally striated, Jacobi ellipsoid, let us define a related, striated sphere and 
focaloidally striated ellipsoid, as matter distributions with same scaled density profile, total 
mass, major equatorial semiaxis and, in the latter case, same boundary. Owing to MacLaurin's 
theorem (e.g., Chandrasekhar 1969, Chap. 3, §19), the gravitational force induced by the stri- 
ated sphere and the focaloidally striated ellipsoid at the end of the major equatorial semiaxis, 
is: 

F G (oi,0,0) = -27rGp 10 (Ai) 10 ai ; (62) 
where p is the mean density, A\ is a shape factor in the range < A\ < 2/3 with the upper 
limit related to spherical configurations, and w — sph, foe, for the striated sphere and the 
focaloidally striated ellipsoid, respectively. 
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The balance between gravitational and centrifugal force at the end of the major equatorial 
semiaxis, reads: 

- 2TrGp w (A 1 ) w a 1 + fi 2 cti = ; (63) 
which, using Eq. ( I37i . may be cast under the equivalent form: 

n 2 

{v eq ) w = n „_ = {Ai) w ; (64) 
ZirLrpw 

where the index, eq, denotes the rotation parameter, v, related to rotational equilibrium at the 
end of major equatorial semiaxis. The above results allow the validity of the relation: 

(Veq)foc < V eq < (v eq ) sph \ (65) 

where v eq is the critical value related to the homeoidally striated, Jacobi ellipsoid. Then v > 
(v e q) foe and v < (v eq ) sp h make a sufficient and a necessary condition, respectively, for the 
occurrence of rotational equilibrium at the end of major equatorial semiaxis, in homeoidally 
striated, Jacobi ellipsoids. 

Owing to MacLaurin's theorem, the condition of rotational equilibrium, Eq. 1641 . for 
focaloidally striated ellipsoids coincides with its counterpart related to homogeneous ellip- 
soids with equal mean density, axis ratios, and velocity field. On the other hand, the rota- 
tion parameter, v, in the special case of focaloidally striated ellipsoids, is also expressed by 
Eq. J40l >. particularized to Jacobi ellipsoids, owing to MacLaurin's theorem. A comparison 
with Eq. i64l shows that rotational equilibrium in focaloidally striated ellipsoids occurs only 
for flat configurations, 631 = 0. Accordingly, Eq. j65t reduces to: 

< v eq < ~ ; (66) 
owing to Eq. ( 1641 and (Ai) sp h = 2/3, as reported in Appendix B. 



4. Transitions from and towards homeoidally striated density profiles 

As an application of the results of Sect. 3, let us take into consideration the transition from an 
initial configuration, with assigned mass, angular momentum, and different kinds of energy, 
towards a final configuration with assigned mass, angular momentum, and different kinds of 
energy tensors. Let both initial and final configurations be homeoidally striated ellipsoids, with 
assigned density profile and rotational velocity profile of the kind discussed in Sect.|2] The 
above mentioned assumption is only for sake of simplicity. It is worth emphasizing that the 
procedure outlined in this Section, well holds for any density profile with triplanar symmetry, 
i.e. symmetric with respect to the principal planes. 

The original method conceived by Thuan & Gott (1975) is aimed to be generalized on 
many respects, namely (i) from spheroidal to ellipsoidal configurations; (ii) from isotropic to 
anisotropic, residual velocity distributions; (iii) from relaxed to unrelaxed, final configurations; 
and, in the limit of axisymmetric configurations (iv) from rigid to differential rotation. 

In the following, first the general case shall be dealt with, in connection with the transitions 
of the kind considered, and then a few interesting special cases will be analysed. 
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4.1. The general case 

The investigation of unrelaxed configurations, for which virial equilibrium does not coincide 
with dynamical (or hydrostatic) equilibrium, goes beyond a purely academic interest. For in- 
stance, (non interacting) galaxies and cluster of galaxies (where galaxies are still falling in 
from outer regions) make examples of (approximately) relaxed and unrelaxed configurations, 
respectively. 

In general, the tensor components, ( pp — ( pp /( and (E res ) pq , which appear in Eq. ( I26K 
have to be considered as effective, instead of actual, anisotropy parameters and residual 
kinetic-energy tensor, respectively. Then the coincidence of virial equilibrium with dynam- 
ical (or hydrostatic) equilibrium is restricted to the special choice £ = 1 which, on the other 
hand, makes a necessary condition only. It is worth recalling that the residual kinetic energy, 
E res , has to be intended as related to all motions different from systematic rotation around a 
fixed axis, such as random motions, streaming flows, coherent oscillations, and so on. With 
regard to the initial configuration, the residual kinetic energy will be expressed as the sum of 
two contributions: one, owing to systematic motions (mainly radial expansion), and one other, 
owing to random motions. 

For sake of brevity, let us define the tensors, together with the related scalar quantities: 

S pq = v se l{B se i) pq ; S — v se iB se i ; (67) 
T^-pq v rarn [B ra7n } p q , 7Z. v raril B rarn , (68) 

and let a prime, and its absence, denote the initial and the final configuration, respectively. The 
effective anisotropy parameters, Q pp , expressed by Eqs. J47> . and the rotation parameters, h, 
£ ro t, v, A, and Xv, expressed by Eqs. J44L ( 145 1 . J46i . J54> . and d57i . respectively, are written 
in compact notation in Appendix C, using Eqs. l|63 and (|68}. 

Owing to Eqs. @, J22> . J67L and J68b . the self-potential and rotational energy related to 
the final and initial configuration, respectively, are: 

E sel = -9^L S ; E' sel = ; (69) 

J 2 J' 2 
Erot = M^ 71 ; K °* = W^ U ' ; (?0) 
and the identity: 

p _ E rot E' rot E' sel , _ 

owing to Eqs. d23l >. ( I69l >. and J70l >. may be cast under the equivalent form: 



0M ^ ^ r 



oi /3j S' Tl S' rc 



(71a) 



Pj = j f ; Pm = ^j ; (7ib) 

involving dimensionless quantities only. 

On the other hand, the identity: 

E , E . , , , 

E = E se i + E, r ot + E res = — E = —(E' sel + E rot + E pec + E osc ) ; 
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by use of Eqs. i24\ . and i26\ . may be cast under the form 

1 E se i 



E se l E se i 
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(72a) 
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where E osc , E pec , represent the kinetic energy of systematic radial and non systematic mo 
tions, respectively. The combination of Eqs. i69\ , ilOi , and i72i yields after some algebra: 

^> £ ^rot — Pe(1 — £ ro t ~ &pec ~ £osc) = i 

involving dimensionless quantities only. 

The combination of Eqs. ( I34L ( I67L J68t . J69t . and ( 17m yields after some algebra 
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which, owing to Eq. ( 145 > . is equivalent to: 
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and Eq. ( I35> may be written under the equivalent form: 

C33 [<5>ii72.22 — <5>2272n] = £33 [C11 72.22 — C22T2.11] ; (76) 

where Eqs. J67l > and d68i have also been used. The substitution of Eq. ill at into I73i produces 
a second-degree equation in £' rot , as: 
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the (reduced) discriminant of this equation is: 
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with regard to a transition from an initial to a final configuration, where all the parameters 
which appear in Eq. (1781 are specified, except the axis ratios, £21 and £31, the condition A = 
represents, via Eq. (17 8> . a curve in the (Oe2i£3i) plane. The transition is forbidden for all 
values of the axis ratios, which make a negative discriminant, and then imaginary solutions. 
The solutions of Eq. (I77a> are: 



£' 



3 



72/ 
72 



2C - 1 C - 1 



2C 



c 



1/2 



(79) 



and the combination of Eqs. (171 at and i79i yields: 
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where the upper and lower sign correspond to the upper and lower sign in Eq. i79\ . It is worth 
noticing that the rotation parameter, £' rot , and the axis ratio, ai/a[, are left unchanged for 
different departures from mass, angular momentum, and energy conservation, provided the 
ratios: 

Ps= w~ E ' Pa= k ] (81) 

do not vary. 

The explicit expression of Eq. ( I76> defines a curve in the (Oe2i£3i) plane, where the axis 
ratios of the final configuration must necessarily lie. The intersection between the curves, 
represented by Eqs. H6\ and d78l >. yields a point where the axis ratios of the limiting, final 
configuration lies, related to a null discriminant. Final configurations with axis ratios, related 
to a negative discriminant, cannot occur. 

4.2. The limiting case c — > 

The limiting case c — > 0, via Eqs. J77l i. corresponds to £' rot = i.e. nonrotating initial config- 
urations, and to £' rot = 2b = 1 — £' osc — £' pec , i.e. unbound configurations. Further inspection 
to Eq. J77cl > discloses that the right hand-side member may tend to zero in a twofold manner, 
through either the ratio £ ro tf ' 0j or the difference within square brakets. 

With regard to the former alternative, the combination of Eqs. J23bt . J71bt . and (1741 yields: 

£rot _ {J') 2 

ft SGM s ai ' ( ' 

which tends to zero for nonrotating initial configurations and/or infinitely extended final con- 
figurations. 

With regard to the latter alternative, the difference within square brakets in Eq. fl77cb tends 

to zero provided the rotation parameter, £ ro t, tends to the value: 

1 1 - 2C 1 
£rot= 2^7; i °<C< g 5 < 83 ) 

where the condition on the virial index, £, is necessary to avoid negative values of the rotation 
parameter in the case under discussion. The combination of Eqs. 1241 . J26b . and 1831 yields 
a second-degree equation involving £ rot as unknown and £ res as parameter, where £ res = 
£ OS c + £pec in the case under discussion. The solutions of the above mentioned equation are: 

£rot ^ j £rot ^ ^res ) (84) 

and the substitution into Eq. (1831 yields: 

1 2£ res 

where the former solution, via Eqs. i24l and ( I26t . corresponds to £ res =3/4, and then to an 
unbound i.e. infinitely extended final configuration, as £ rot + £ res = 1. On the other hand, the 
latter solution clearly represents a relaxed configuration, as £ rot + £ res =0.5, which would 
necessarily imply ( — 1, contrary to what has been found. Then the limiting case c — > is 
attained by making the ratio £ rot / 0j tend to zero. 
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4.3. Virialized, final configurations 

In the special case of relaxed, final configurations, £ = 1, virial equilibrium coincides with 
dynamical (or hydrostatic) equilibrium. Accordingly, Eq. (I77c> translates into: 
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similarly, Eqs. i!9i and d80i reduce to: 
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where the upper and lower sign must correspond in both equations. 



4.4. Nonrotating, initial configurations 

In the special case of initial configurations with vanishing angular momentum, J' — > 0, the 
parameter, PJ , has necessarily to remain finite. Accordingly, the combination of Eqs. ( I23bt . 
<E3B, and CD yields: 
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where the product on the left-hand side of Eq. ( I90> is infinitesimal of higher order with respect 
to its counterpart related to Eq. ( I89> . Then Eq. (I77at reduces to a first-degree equation, and 
Eqs. |79) and (|80j reduce to: 
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(92) 



respectively. It is worth noticing that unbound configurations correspond to b = 0, via 
Eqs. J72bl > and dTTbt . 

If, in addition, the rotation energy of the final configuration is negligible in respect of the 
self -potential energy, £ ro t — > 0, Eqs. ( I9U and d92i reduce to: 



Pj rot ~ 2(3 E {S'J K rot 
ai _ 1 P 2 M S 2C - 1 1 
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(93) 
(94) 



where both sides of Eq. (1931 have necessarily to be infinitesimal of the same order. 

In the special case of a relaxed, final configuration, which necessarily implies £ = 1 
according to Eqs. \2A\ and J26i . Eqs. J87t and J88I reduce to Eqs. \9\\ and ( I92l i. respectively, 



(95) 
(96) 



particularized to the special case, ( = 1. The result is: 
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where the trend is the same as shown above. 
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4.5. Application to dark matter haloes 

Dark matter haloes make a natural application of the procedure outlined above, as they are the 
counterpart of elliptical galaxies passing from baryonic to dark matter universes. Aiming to see 
nothing but the method at work, our attention shall be restricted here to (i) CDM cosmological 
models and (ii) (non baryonic) dark matter, which is equivalent to assume a homeoidally 
striated density profile also for baryonic matter where, in addition, the isopycnic surfaces are 
similar and similarly placed with respect to the boundary (e.g., CM03). 

The current procedure has a main advantage on numerical simulations, that the initial 
configuration may be taken as an overdensity region at recombination epoch. In this case, 
the dominant energies are expected to be the potential and the kinetic one, the latter related 
to systematic, expansion motions, and the mass distribution may safely be thought of as a 
homogeneous sphere. Computations usually start much later, owing to their high cost and time 
machine. The choice of recombination epoch as starting time makes the initial configuration 
virtually independent of the cosmological model, at least for sufficiently large, present-day 
density parameters (f2o > 0.1), where the related universe may be considered flat to a good 
extent i.e. |1 - f2~ e * | < 1 (e.g., Zeldovich & Novikov 1982, Chap. Ill, §4; Caimmi & Marmo 
2004). 

Let us represent overdense regions as made of two subsystems with coinciding volume, 
namely (i) the unperturbed, cosmological fluid, with density equal to the mean density of 
the universe, ph, and (ii) the density perturbation, with mean density equal to the difference 
between the overdensity and the mean density of the universe, bp — p — pt- At recombination 
epoch, overdense regions may safely be thought of as homogeneous and spherical. This is why 
the overdensity index, S rec = Sp/p^, cannot exceed (in absolute value) a few percents at that 
time (e.g., Gunn 1987). In addition, tidal distorsions from sphericity are owing to neighbouring 
density perturbations, while tidal distorsions towards sphericity are owing to the unperturbed, 
cosmological fluid. 

In general, the potential-energy tensor of an overdensity region is the sum of three contri- 
butions, namely: (i) the self-energy tensor of the overdense region, (ii) the tidal-energy tensor 
induced by the unperturbed, cosmological fluid outside the overdense region, and (iii) the tidal- 
energy tensor induced by neighbouring density perturbations. For reasons mentioned above, 
the last two contributions may safely be neglected in the case under discussion. Accordingly, 
the combination of Eqs. @, Q, and J67> : J22l > and d68i : yields: 



respectively, as v se i =3/10 and A\ = A2 — A3 = 2/3 in the case under discussion of a ho- 
mogeneous sphere; and r\ anm — 2r) ro t — 1, v anm = v rot = Vi nr = 1/5 for a homogeneous, 
rigidly rotating sphere. 



S'pq ~ ^ S pq 1 S' ~ - ; 
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^'pg = gSpqO- - spa) ; ft — - ; 
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The kinetic-energy tensor of systematic, expansion motions, is: 

(Ksc) Pq = ^S pq M'a' 2 ; E' osc = ^M'a' 2 ; (99) 

in the case under discussion of a homogeneous sphere. The further assumption that, at recom- 
bination epoch, the overdense region expands at the same rate as the unperturbed, cosmologi- 
cal fluid, implies the validity of the relation: 

& = (a,h)rec = H rec {a,h)rec = H rec a' ; (100) 

where ah is the scale factor of the universe, H is the Hubble parameter, and the index rec 
denotes the recombination epoch. 

At sufficiently early times i.e. high redshifts, in particular at the recombination epoch 
(z rec ~ 1400), the following relations hold (e.g., Zeldovich & Novikov 1982, Chap. Ill, §4): 

n- 1 (z) = l+ 1 1 n "° ; Ho > 0.1; (101) 

Phiz) = + ; |i - « i ; d02) 
h(z) = Hi = h (i + zf' 2 ■ |i - rr 1 ! « 1 ; (103) 

where ft is the density parameter (to be not confused with the angular velocity), z the red- 
shift, t the time, and the index denotes the present-day epoch. The definition of the density 
parameter reads: 

n = r-r— ; toW = jr? > ( 104 ) 

where (ph)crit is the critical density i.e. the density of the Hubble flow in an Einstein-de Sitter 
universe. 

By use of the definition of mean density inside an overdense region, p = + Sp = 
p h (l +6) = 3M h (l + d)/(4nal), and Eqs. (fTOOt . ST021 . and (fT03l . the following relations 
are derived after some algebra: 
2GM' 



H 2 ^ 



1/3 1 

1 -\- z re 



(105) 



i 1/3 

a! = H (l + z rec fl 2 o! = ( 2 JfM2k\ (i + Zrec) i/2 . (106) 



V 1 +S r 

and the combination of Eqs. (I99i and JlOOi reads: 

(Ksc) Pq = ^5 pq H 2 rec M'a' 2 ; E' osc = ^H 2 rec M'a' 2 ; (107) 

in the case under consideration of a homogeneous sphere expanding at the same rate as the 
universe. The combination of Eqs. J67i . ( I69> . (I72b> . (I104t and dl07> yields: 

Cc= [Orec(l + ^ e c)] _1 ; (108) 

which has necessarily to be less than unity for bound density perturbations. 

Accordingly, the sum of initial systematic rotational and peculiar energy, E' rot + E' pec , 
cannot exceed ~E' pot — E' osc for bound configurations, i.e. £' rot + £' pec < 1 — £' osc . Keeping 
in mind that il rec may safely be put equal to unity and the overdensity index, S rec , amounts 
to a few percent, it can be seen from Eq. dl08t that the sum of parameters, £' rot + £' pec , cannot 
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exceed a few percent, or in other words the sum of initial rotational and peculiar energy is two 
orders of magnitude lower than the initial expansion energy. Then the knowledge of the initial 
rotational and peculiar velocity field is not essential, and for this reason an initial rigid rotation 
is assumed, together with a null, initial, peculiar velocity field. 

With regard to the final configuration, the simplest case of mass conservation, M = M' 
or (3m — 1, and energy conservation, E = E' or (3e = 1, makes a valid approximation for 
isolated, overdense regions. On the other hand, merger accretion necessarily implies (3m > 1 
and (3e > 1. In addition, (complete) virialization, E se i + 2(E rot + E res ) — 0, may safely 
be assumed in dealing with density perturbations of galactic size, while it must be released 
for larger dimensions i.e. clusters of galaxies, where the inner region is relaxed but the outer 
shells are still infalling. 

For representing density profiles of virialized, dark matter haloes, expressed by Eqs. Q, 
two alternatives will be exploited. One, first proposed by Navarro et al. (1995; 1996; 1997), 
reads: 

f(0 = ; (109) 

and will be quoted in the following as the NFW density profile. One other, first proposed by 
Moore et al. (1998; 1999), reads: 

^) = g 3/ 2(1 2 +C 3/ 2) i ( 110 ) 

and will be quoted in the following as the MOA density profile. In both cases, the truncated 
scaled radius is determined by fitting the results of high-resolution simulations in a CDM cos- 
mological model (Fukushige & Makino 2001) to the selected density profile. For more details, 
see CM03; Caimmi & Marmo (2004). In addition, both rigid rotation (triaxial configurations) 
and constant rotational velocity on the equatorial plane (axisymmetric configurations) are as- 
sumed. In the latter alternative, both rigid (r) anm — 2rj rot — 1) and differential (to ensure 
constant rotational velocity everywhere: r) anm — 3it/8, r\ rot = 3/4) rotation of the generic, 
isopycnic surface, are considered. The values of the truncated, scaled radius, and some profile 
parameters, expressed in Sect. [2] and particularized to both NFW and MOA density profiles, 
are listed in Tab.[5] 

At this stage, the calculations are allowed to start, following the procedure outlined below. 

(i) Define the values of the input parameters: the axis ratios, e' 21 and e' 31 , the profile parame- 
ters, v' sel and v' ram , the overdensity index, <5 rec , the degree of departure from mass conser- 
vation, (3m, energy conservation, (3e, angular momentum conservation, (3j, and the virial 
index, £, the last in connection with the final configuration. 

(ii) Select a density profile for the final configuration, and calculate the profile parameters, 

^masi ^inr-> ^sel- 

(iii) Select a velocity profile (of the kind considered in Sect.|2ji for the final configuration, 
and calculate the profile parameters, v anm , v rot , v ram , and the shape parameters, rj anm , 

Vrot- 
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(iv) Select a sequence of homeoidally striated Jacobi ellipsoids (which represent final con- 
figurations), where the generalized and effective anisotropy parameters, ( p and ( pp , are 
kept constant. 

(v) Select a value of the axis ratio, £21, of a configuration on the sequence under considera- 
tion. 

(vi) Calculate the value of the remaining axis ratio, £31. 

(vii) Calculate the value of the rotation parameters, h, £ ro t, v, A, and Xv> defined by 
Eqs. C3, (Q, iC4t . dC5t . and Ol . respectively. 

(viii) Calculate the value of the ratios, £' rot / (1 — £' osc ), and ai/a[, defined by Eqs. J77i . (179b . 
{80j, and JTOSl . 

(ix) Select another configuration on the sequence under consideration, and return to (v), 
unless the whole sequence is covered. 

(x) Select another sequence of homeoidally striated Jacobi ellipsoids, where the generalized 
and effective anisotropy parameters, ( pp and (pp, are kept constant, and return to (v). 

(xi) Select another velocity profile for the final configuration, and return to (iv), unless the 
range of interest is covered. 

(xii) Select another density profile for the final configuration, and return to (iii), unless the 
range of interest is covered. 

(xiv) Select another value of the virial index, £ (in connection with the final configuration), 
and return to (ii), unless the range of interest is covered. 

4.6. Results and discussion 

Our attention shall be limited to dark matter haloes hosting giant galaxies, with masses of 
about 10 12 m Q . Accordingly, we assume S rec — 0.015 as a typical overdensity index of the 
initial configuration, and a virialized final configuration i.e. £ = 1. In addition, the following 
range of parameters, (3g and j3 a , is considered: 

— > B e > — !— ; 100 < 8 a < 14400 ; (111) 
100 ~ H ~ 14400 v > 

which includes the case /3m = 0e = 1> allowing (3j lie within the range: 

10 < (3 j < 120 ; (112) 



owing to Eq. (18 II . For sake of simplicity, in what follows the last choice of (/3m, Pe, will 
be adopted, but keeping in mind that it is equivalent to any other choice satisfying Eqs. (II 1 li . 
for what will be analysed and discussed. 

In particular, three sequences of density profiles with constant anisotropy parameters, £ pp , 
are considered, starting from a nonrotating, axisymmetric configuration equal to an oblate 
(eai = 0.5; C11 = 0.3910; C33 - 0.2180), round (e 3 i = 1; C11 = C33 = 1/3), and prolate 
(e 3 i = 2; Cn = 0.2717; C 33 = 0.45 66) spheroid, respectively. 

As in Thuan & Gott (1975) original approach, a limiting configuration exists, which makes 
the discriminant, expressed by Eq. (I78l >. equal to zero. More flattened configurations are not 
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Fig. 5. The axis ratios of the limiting configuration as a function of the parameter, f3j, for 
sequences of virialized density profiles with constant anisotropy parameters, starting from 
axisymmetric, nonrotating configurations, defined as: O (£31 = 0.5; (11 = 0.3910; £33 = 
0.2180); S (e 3 i = 1; C11 = C33 = 1/3); P (e 3 i = 2; (n = 0.2717; C33 = 0.4566). The 
polar and the equatorial axis ratio, £31 and £21, are represented by full and dashed curves, re- 
spectively. Configurations for which bifurcation from axisymmetric to triaxial density profiles 
occurs, are marked by filled circles. Pentagon skeletons mark configurations for which the va- 
lidity of a necessary condition for the occurrence of centrifugal support at the ends of the major 
axis, is first satisfied. The density profile and the kind of rotation are indicated in each panel. 
Configurations related to constant rotational velocity on the equatorial plane, v(ir/2) = const, 
and rigid rotation of isopycnic surfaces, have no physical counterpart for triaxial boundaries. 

allowed for the related sequence of density profiles and choice of (/3m, Pe, A/)- The depen- 
dence of the axis ratios of the limiting configuration on the parameter, f3j, related to the three 
sequences mentioned above, is represented in Fig.[5]for NFW and MOA virialized density 
profiles, characterized by either rigid rotation or constant rotational velocity on the equato- 
rial plane and rigid rotation of isopycnic surfaces. Configurations for which bifurcation from 
axisymmetric to triaxial density profiles occurs, are marked by filled circles. Pentagon skele- 
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tons mark configurations for which the validity of a necessary condition for the occurrence of 
centrifugal support at the ends of the major axis, according to Eq. (1651 . is first satisfied. 

As low values of (3j correspond to bounder final configurations and vice versa, the axis 
ratio of the limiting configurations are strongly dependent on the sequence under consideration 
for low values of j3j and vice versa. In other words, a large amount of rotational energy erases 
the memory of the initial shape, different for different sequences. A complete dissipation of 
angular momentum, (3j — 0, yields nonrotating final configurations. In addition, it can be seen 
from Fig. |5] that a conservation of angular momentum, (3j = 1, implies final configurations 
lying near the nonrotating limit. 

The density profile and the kind of rotation are indicated in each panel of Fig. [5] Configura- 
tions related to constant rotational velocity on the equatorial plane, v(tt/2) = const, and rigid 
rotation of isopycnic surfaces, have no physical counterpart for triaxial boundaries, but they 
may be of mathematical interest for comparison with solid-body rotating mass distributions 
with equal density profile. On the other hand, a change in density profile affects the curves to 
a smaller extent for configurations with constant rotational velocity on the equatorial plane, 
with respect to rigidly rotating configurations, as expected. 

From this point on, our attention shall be focused on two special values of (3j, namely: 

13 j = 60 ; 0.20 £ (e 31 ) Um £ 0.40 ; 0.30 £ {e 2X ) Um £ 0.50 ; (113) 
0j = 120 ; 0.05 £ (e 31 ) iim ~ 0.16 ; 0.06 £ (e 21 ) Hm £ 0.15 ; (114) 

in connection with NFW density profiles in rigid rotation. An inspection to Fig.[5]shows that 
the validity of a necessary condition for the occurrence of centrifugal support at the ends of 
the major axis, according to Eq. i65\ . may be violated. 

The axis ratios of the final configuration, £31 and 621 (upper panels) and the axis ratio 
of final to initial configuration, a\ja! x (lower panels), as a function of the parameter, \£ — 
£' rot /(l — £' osc ), are plotted in Fig.[6]for three sequences of density profiles, O, S, and P, 
considered earlier, and /3j = 60 (left panels) and (3j — 120 (right panels). Each curve on the 
upper panels is symmetric with respect to a vertical axis, \£ — 0-5, and is characterized by 
three extremal points, two maxima and one minimum. The maxima are related to the starting 
configuration of the sequence, with finite and infinite axes, respectively. On the other hand, 
the minimum corresponds to the limiting configuration, which exhibits the lowest values of 
the axis ratios, allowed for the sequence under consideration. Flat configurations occur only 
in the limit j3j — > +00 unless centrifugal support takes place at the ends of the major axis, 
from which two flows of dark matter freely stream (see e.g., Jeans 1929, Chap. XIII, § 301, 
with regard to gaseous nebulae). In addition, larger values of f3j make different curves closer 
each to the other, and vice versa. The limit, \£ ~ > 1> corresponds to unbound configurations 
(E = 0) which extend to infinite. 

The rotation parameters, £ ro t, h, A, and v, (from top left in clockwise sense), and \v, 
related to the final configuration, as a function of the axis ratios, £31 and £21, are plotted in 
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Fig. 6. The axis ratios of the final configuration, 631 and £21 (upper panels) and the decimal 
logarithm of the axis ratio, ai/a[, of final to initial configuration (lower panels), as a function 
of the parameter, \£ = £rot/ (1 — £'osc)> f° r sequences of virialized density profiles, O, S, and 
P, considered in Fig.|5] top left panel, related to j3j = 60 (left panels) and /3j — 120 (right 
panels). The curves on the upper panels are symmetric with respect to a vertical axis, \£ — 0-5. 
The point of minimum of each curve of the upper panels corresponds to the related, limiting 
value of the axis ratio, i.e. the lower value allowed along the sequence under consideration. 
Other captions as in Fig.|5] The limit, xs —> 1> is related to unbound configurations (E = 0) 
which extend to infinite. 

Figs. and [8] respectively, for three sequences of density profiles, O, S, and P, considered 
earlier. The curves are independent of the parameter, (3j, but the limiting values of the axis 
ratios, related to the selected value of (3j, cannot be exceeded moving from the right to the left 
along each curve. The horizontal lines correspond to a mean value, A = 0.0421, and related 
variations within one rms error, A + = 0.133 and A~ = 0.0133, consistent with the results of 
high-resolution numerical simulations (e.g., Bullock et al. 2001; Vitvitska et al. 2002). 

Leaving aside peculiar situations such as the occurrence of a major merger near the end 
of evolution, the main results of high-resolution simulations (e.g., Fukushige & Makino 2001; 
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Fig. 7. The rotation parameters, £ ro t, h, A, and v, (from top left in clockwise sense), related 
to the final configuration, as a function of the axis ratios, 631 (full curves) and £21 (dashed 
curves), for sequences of density profiles, O, S, and P, considered in Fig.|5] top left panel. 
The curves are independent of the parameter, 0j, but the limiting values of the axis ratios, 
related to (3j, cannot be exceeded moving from the right to the left along each curve. The 
horizontal lines correspond to a mean value, A = 0.0421, and related variations within one rms 
error, A + = 0.133 and A~ = 0.0133, consistent with results from high-resolution numerical 
simulations. Other captions as in Fig.[5] 

Klypin et al. 2001; Bullock et al. 2001; Vitvitska et al. 2002) may be summarized as follows, 
to a good extent: (i) a substantial fraction of the initial density perturbation gets virialized; (ii) 
the density profile may be considered as self-similar and universal; and (iii) the distribution 
of the rotation parameter, A, is lognormal with mean value, 0.03 ~ A ~ 0.05, and rms error, 
0.5 ~ cri og x ~ 0.7, in particular A = 0.0421, cr log A = 0.500. 

Accordingly, simulated density profiles with assigned mass, boundary, and rotation param- 
eter, may be compared with their counterparts deduced from the current model. In particular, 
an inspection of Figs.[6]and0shows that, in the range 0.013 < A < 0.133, dark matter haloes 
are close to the starting point of the selected sequence i.e. the nonrotating configuration. Then 
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Fig. 8. The rotation parameter, %„, related to the final configuration, as a function of the axis 
ratios, £31 (upper panel) and £21 (lower panel), for sequences of density profiles, O, S, and P, 
considered in Fig.[5] top left panel. The curves are independent of the parameter, 0j, but the 
limiting values of the axis ratios, related to j3j, cannot be exceeded moving from the right to 
the left along each curve. Other captions as in Fig.[5] 

the shape of dark matter haloes is decided mainly by the amount of anisotropy in residual ve- 
locity distribution, while the contribution of rotation has only a minor effect on the meridional 
plane, and no effect on the equatorial plane. This is why bifurcation points occur for larger 
values of A. In other words, dark matter haloes may exhibit any shape, provided it is owing 
mainly to anisotropy in residual velocity distribution and only weakly to systematic rotation. 

The analysis of elongated configurations which rotate around the major axis goes beyond 
a purely academic interest. Empirical evidence in favour of the existence of galaxies with 
prolate stellar structure, cut equatorially by a gaseous lane, was found long time ago (Bertola 
& Galletta 1978). The hosting dark haloes could also be prolate and spinning around the 
symmetry axis. A recent investigation on polar ring galaxies has indicated that the hosting 
dark haloes are most likely elongated towards the polar ring plane (Iodice et al. 2003), where 
the rotation axis of the inner body lies. 
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5. Conclusion 

The formulation of the tensor virial equations has been generalized to unrelaxed configura- 
tions, where virial equilibrium does not coincide with dynamical (or hydrostatic) equilibrium. 
Special classes of homeoidally striated ellipsoids have been shown to exhibit similar properties 
as Jacobi ellipsoids, and to reduce to Jacobi ellipsoids in the limiting situation of homogeneous 
matter distribution, rigid rotation, and isotropic residual velocity distribution. Accordingly, the 
above mentioned density profiles have been defined as homeoidally striated, Jacobi ellipsoids. 
Further investigation has been devoted to the generation of sequences of virial equilibrium 
configurations, and sequences with constant rotation and anisotropy parameters have been 
studied with more detail, including both flattened and elongated triaxial (in particular, both 
oblate and prolate axisymmetric) configurations, and the determination of the related bifurca- 
tion points. The explicit expression of a number of rotation parameters, used in literature, and 
the ratio of rms rotation to rms residual velocity, have also been calculated. 

Thuan & Gott (1975) procedure has been generalized to homeoidally striated, Jacobi el- 
lipsoids, and a few limiting situations have also been analysed. An application has been made 
to dark matter haloes, and calculations have been performed for masses M w lO 12 m i.e. 
haloes hosting giant galaxies. The results have been shown to be qualitatively similar to what 
was found by Thuan & Gott (1975), in particular concerning the occurrence of limiting axis 
ratios, below which no configuration is allowed for the sequence under consideration. The 
dependence of the limiting axis ratios on two parameters, related to the degree of change in 
mass, total energy, and angular momentum, has been illustrated in connection with NFW and 
MOA density profiles, rotating with constant either angular or linear (on the equatorial plane) 
velocity, for three special sequences starting from an oblate (631 = 0.5), round (£31 = 1), and 
prolate (£31 = 2), nonrotating configuration, respectively. 

The dependence of the axis ratios, £31 and £21, on the parameter, \£ = £' ro t/ (1 — ^osc)' 
has been analysed in connection with rigidly rotating, NFW density profiles, for the three 
special sequences mentioned above. The same has been done for the rotation parameters. 
Within the range of the rotation parameter, exp 10 (log A=F(7i og A) = exp 10 (— 1.3761 =f0.500), 
consistent with results from high-resolution numerical simulations, it has been shown that the 
shape of dark matter haloes is mainly decided by the amount of anisotropy in residual velocity 
distribution. On the other hand, the contribution of rotation has exhibited only a minor effect 
on the meridional plane, and no effect on the equatorial plane, as bifurcation points occur 
for larger values of A. To this respect, dark matter haloes have been found to resemble giant 
elliptical galaxies (e.g., Binney 1976; Binney & Tremaine 1987, Chap. 4, § 3). 
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Appendix A: Different normalizations for homeoidally striated density 
profiles 

With regard to homeoidally striated density profiles, defined by Eqs. Q, the mass, the inertia 
tensor, and the self-energy tensor, may be expressed as (CM03): 

M = vl as M^ ; (Al) 

i Pq = <wL Mt (4) 2 ; < A2 ) 

+ G(M^) 2 

(E sel ) pq = -u\ el y 1 (B sel ) pq ; (A3) 
a[ 

where the profile parameters, v\ take the explicit form (CM03): 



vl = 

mas 



~J~F(0d£ ; (A4) 

v\nr = l f 'me^ ; (A5) 



2 7o 



-:.6./ () F (e)de; (A6) 



sel 16 Jo 

and the function, has the definition (Roberts 1962): 

F(0 = 2 jT /(£')£' ; (A7) 
where the function, /(^), corresponds to the density profile, according to Eqs. ([0. It can easily 



be seen that the following relations hold: 

F(E) = ; ^ = -2£/(0 ; (A8) 
and an integration by parts shows that: 

77 — 1 r s 

f(m n ^ = -^r- / F (0C~ 2 ^ ; n > 1 ; (A9) 



o 2 j 

which allows the calculation of the profile parameters. The scaling mass, M\ has the same 
definition in CM03 as its counterpart, Mq, in the current paper, according to Eq.(|6}. Then the 
following relation holds: 

A/t = Mq ■ (A10 ) 
for further details, see CM03. 
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The mass, the inertia tensor, and the self-energy tensor, in the current paper are expressed 
by Eqs.0, 0, and ®, respectively. The comparison with their counterparts, expressed by 
Eqs. (ED, (E3, and dA3l . yields: 



Vrnas — v mas i (All) 
-J ^ ^ • 

M a, 



i (A 12) 



p 

"•«' = "-(^) 3 t ; ,A13) 

where the last two relations, owing to Eqs. dlbl . (0, ( IA10K and iAl 1> . may be written under 

the equivalent form: 

,t 



L/ raas^- i 

vis 



(A14) 



Vsel = 7T^-2 5 (A15) 

[yraas ) 

and the profile parameters, v\ nas , v\ , v* are expressible as functions of the scaled radius, 
3 (CM03). 

The angular-momentum vector and the rotational-energy tensor, for the special rotational 
velocity fields discussed in Sect.|2] may be expressed as (CM03): 

Js = $s3vlnmvlnm Mia t( 1 + e lp)( v tot)p 5 P + 1 + * 5 (A 16) 



{Erot) Pq = 6pq(l - 5 p3 )rit ot vl ot M'< (vl ot ) p 



(A17) 



and the related module (angular momentum) and trace (rotational energy), respectively, as: 
J = 4nm4nmM^al(l + el 1 )(vt ot ) 1 ; (A18) 
Erot = vL4otM*{l + 4i)[(vl at ) 1 ] 2 ; (A19) 
where {v\, ot ) p is the rotational velocity at the top axis, a|, of the reference isopycnic surface. 
The profile parameters, v\ and the shape parameters, rf , take the explicit form (CM03): 

^ 2(C ' 0) /(£)£ 4 d£ ; (A21) 



o O 2 (l,0) 

3 4 /" +7r/2 0(1,61) sin 3 (9d(9 



• (A22) 

4^ ./.^ 0(1,0) (cos 2 9 + e 2 sin 2 (9) 5 / 2 ' 

t 3 4 f+-/ 2 Q 2 (l,6Q Bin 3 gdg . 

™* 8 7_ w/2 O 2 (l,O)(cos 2 + e 2 sin 2 0)5/ 2 ' k } 

where e is the ratio of polar to equatorial axis and 8 is the azimuthal angle. It is worth of 

note that the above expressions for the shape parameters, if, hold for axisymmetric configura- 
tions. The shape parameters, r\, are also expressed by Eqs. JA22t and JA23> . provided 0(1, 9) 
and O(l,0) therein are replaced by 0(3,6) and 0(3, 0), respectively. The extension of the 
above definitions and results to (necessarily) solid-body rotating, triaxial configurations, yields 
Vanm = 1> Vrot = 1/2, according to Eqs. (I7l-(I10||. For further details, see CM03. 

The angular-momentum vector, the rotational-energy tensor, the angular momentum, and 
the rotational energy, in the current paper are expressed by Eqs. Jl 8i . dl9l >. and i20\ . 12 H . 
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respectively. The comparison with their counterparts, Eqs. ( IA16I . 1A17I . and (IA18> . dA19i . 

respectively, yields: 

v _„ t "Lm ^ 4 (vj ot ) P _ „_ 12 . (A2 4) 

ijanm ®"p \Prot )p 

iVot = ^ ot ttt ; P = 1.2 ; (A25) 

where, owing to Eq. ?7 onm = r\\ nm and ?7 rot = ?j| ot . Then Eqs. (IA24> and ( IA25> . by use 
of (llb> . (|2j, and JA10> . take the equivalent form: 

P=l,2; (A26) 



"L (4ot)p _ "L o(i, o) 



l4u,«S («rat)p ^„a S S 2 0(3, 0) 
4* frjjg jjg » 2 (1,0) 1 - 

iVot = -l — t = -: — — 2 ,„ , ; p= 1,2 ; (A27) 

where {v\. ot ) p / {v rot ) p = 1/Hif the system is in rigid rotation, and (v\ ot ) p / '(v rot ) p = 1 if the 
rotational velocity is constant on the equatorial plane. The profile parameters, v\ nm , v\ ot , are 
expressible as functions of the scaled radius, 5 (CM03). 

The above results allow the calculation of the profile parameters, v, appearing in the text. 
The corresponding, explicit expressions, may be seen in Caimmi & Marmo (2004). 



Appendix B: Limiting configurations 

The general expressions of the shape factors, A p , and related quantities, for limiting configu- 
rations of ellipsoidal shape reduce to undetermined expressions of the kind 0/0, oo/ oo, et sim- 
ilia (e.g., Caimmi 1991, 1992, 1995). Then explicit calculations are needed for round, oblate, 
prolate, flat, and oblong configurations, which can be found in the above quoted references. 
Though flat and oblong configurations are dynamically unstable (e.g., Jeans 1929, Chap. IX), 
they shall also be included for sake of completeness. Accordingly, the expression of the shape 
factors in the self-energy tensor, (B se i) pp — e P 2€ P 3,A p , see Eq. (JSJl, must be determined even 
for flat and oblong configurations. 

The shape factors, A p , depend on the axis ratios, via elliptic integrals of first and second 
kind (e.g., Caimmi 1992): 

F(0,p) = 



o yl — p 2 sin 



2 



E((3,p) = [ Jl-p 2 sin 2 ^ ; (Bla) 
Jo 



e-21 



13 = arccose 3 i ; p = — ; e pq = Jl - e 2 ; (Bib) 

e3i v 

where (3 = n/2 for both flat and oblong configurations, while p = e%\ and p — 1 for flat and 
oblong configurations, respectively. 

With regard to the self-energy tensor, by use of the general expression of the shape factors, 
A p , derived in Caimmi (1992; therein defined as a p ), the flat limit (£21 > 0, £31 — > 0) reads: 

1 "' ' ' = 61 ; (B2a) 

= h ; (B2b) 
; (B2c) 



(B se i)n = ei 2 ei 3 Ai = F (J^, e 2 %j - E (J^, e 

{B se l)22 = £22623^2 = E \ 2 ' e21 J _ £ 21 i? ( 2 ' 



2 

{B S el)d,3 = £32633^3 



e2i H e 2 i 



21 L KZ ' £ 21 
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and the special case of the oblate, flat limit (£21 = 1, £31 — * 0) reads: 



(B ae i)u = [B se i) 2 2 = lim (e D V) = - 



(B se i) 33 = lim. (e„7) = 



£o — £31 — £32 ; 



a 



lim = lim 

£21— >l e 2 i— >l 



7 



lim A3 

621— >1 



(B3a) 
(B3b) 



for further details see e.g., Caimmi (1991). 

Oblong configurations may be conceived in a twofold manner, according if £21 — > £31 = 
(flat-oblong, hereafter quoted as f-oblong) or £21 = £31 — * (prolate-oblong, hereafter quoted 
as p-oblong). In any case, f3 — ► w/2, p — ► 1, e2i — * 1, and Eqs. JB 1 at reduce to: 

1 + tan x 



lim f(-, 

■21^0 V2 



f J 21 



lim In ■ 



(I 



1 



(B4) 



x^7r/4 1 — tana; 

in addition, using the series development of the complete elliptic integral of first kind and 
Wallis' product (e.g., Spiegel 1968, §§ 34.2 and 38.9) yields: 



2 nfir \ 7T .. [l-3-...-(2n-l) 
lim n 4 1 F[-,e 2 i )=- lim — — i— '- 

with regard to f-oblong configurations, Eqs. (IB2t reduce to: 

{B se i)n = +00 ; (B sei )22 = 2 ; (S sei ) 33 = : 



= lim 



-00 2n + 1 



= 



(B5) 



(B6) 



on the other hand, with regard to p-oblong configurations, the shape factors in the self-energy 
tensor take the expression: 

2 



(B sel )u = lim (e_ 2 7)= lim e 2 

e„-+0 r e— >+oo 



{B se i) 22 = {B se i) 33 = lima = lim 



1 



arcsinh \/£ 2 — 1 
£ f= — 1 



1 



arcsinh y/e 2 — 1 



= 1 



(B7a) 



(B7b) 



£31 — £21 ; 



lim Ai ; 

621— *€31 



7 = lim A 2 = lim A 3 ; (B7c) 

621— *€3i 621—^631 



for further details see Caimmi (1991). 

The comparison between Eqs. dB6l and (IB7> discloses that the shape factors in the self- 
energy tensor exhibit a discontinuity passing from f-oblong to p-oblong configurations (e.g., 
Caimmi 1993b). The values of the axis ratios, £21 and £31, the shape factors, A p and (B se i) pp , 
related to limiting configurations of ellipsoidal shape, are listed in Tab. lB II where the upper 
part, containing the shape factors, A p , has been taken from Caimmi (1995). The range of 
validity of the independent variables, £/, e p , e , and e, and the definition of the functions, a, 
7, are: 



< £/ < 1 ; < e p < 1 ; < £ G < 1 : 
< £ < 1 , oblate; 1 < e < +00 , prolate; 



(B8a) 
(B8b) 




2U/2 



arcsin(l — e ) 

arcsinh(£ 2 — 1) 

" (£ 2 - 1) 1/2 



1/2 



oblate; 
; prolate; 



(B9a) 



R. Caimmi and C. Marmo: homeoidally striated density profiles 



261 



7 



1 -e z 
2 



1-e 



arcsin(l — e ) 



2U/2 



(l-e 



2U72" 



arcsinh(e 2 — 1 



,1/2 



- 1 



oblate; 
prolate; 



(B9b) 



(e 2 - if^ 

and the functions, b\ and 62, are expressed by Eqs. JB2a> and (IB2bt . For further details, see 
Caimmi (1995). 



Appendix C: Expression of anisotropy and rotation parameters in 
compact notation 

The effective anisotropy parameters, ( pp , expressed by Eqs. ( I47K and the rotation parameters, 
h, £ rot , v, A, and \v, expressed by Eqs. ( l44l . J45l l. j46l >. d54l . and ( BTt . respectively, may be 
written in compact notation, using Eqs. J67l > and J68b . The result is: 



5, 



c 



pp 
5 



7? 

1-2/i^ 








[—!] 









p= 1,2,3 



with regard to anisotropy parameters; 



2ft 



^ C ^33 

C33 <5 



&rot — § 

1 1 

v = — - 



£21^31 



3 rj ro tV ro t 1 + C 2 1 



5 



, C ^33 

C33 5 



1 - 



C £33 
C33 5 

i-C C S33 
C C33 s 



xl = C 



1 



4 ?e 

C33 s 
C £33 

with regard to rotation parameters. 

Finally, the combination of Eqs. dC2l i. dC4l >. jC5> . and jC6l > yields: 

2 1 £2l£31 rry 1 

v = - 5-ft/i ; 

3 r) ro tv rot l + fi 



(CI) 

(C2) 

(C3) 
(C4) 

(C5) 

(C6) 



A = -Sh 



1 - 



21 

1-C C ^33 
C C33 S 



xl = C- 



2Kh 



(C7) 
(C8) 
(C9) 



' S - 2TZh 

which allows a connection between different rotation parameters. 

An inspection of Eq. (ICU shows that the rotation parameter, h, attains the special values 
1 S, 



2 l^qq 
1 S -Sr, 



9 = 1,2 



2ft -ft 



Cpp — i 



p = l,2,3 



/'/< 



h = 



5, 



pp 
5 



p= 1,2,3 ; 



(ClOa) 
(CI Ob) 
(ClOc) 



where the last relation is equivalent to Eq. (|48 
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£31 £21 Srot tlN Vn Ajv XN 



1.000000 


1.00 


0.0000000 


0.0000000 


0.00000000 


0.000000 


0.000000 


0.950000 


1.00 


0.0136423 


0.0554973 


0.00878707 


0.237573 


0.167481 


0.900000 


1.00 


0.0279335 


0.115614 


0.0173421 


0.345874 


0.243255 


0.850000 


1.00 


0.0429239 


0.180831 


0.0256177 


0.436408 


0.306447 


0.800000 


1.00 


0.0586689 


0.251690 


0.0335587 


0.519556 


0.364604 


0.750000 


1.00 


0.0752304 


0.328809 


0.0411011 


0.599400 


0.420843 


0.700000 


1.00 


0.0926773 


0.412889 


0.0481704 


0.678136 


0.476999 


0.650000 


1.00 


0.111086 


0.504734 


0.0546795 


0.757185 


0.534446 


0.600000 


1.00 


0.130544 


0.605265 


0.0605265 


0.837600 


0.594425 


0.582724 


1.00 


0.137528 


0.642207 


0.0623716 


0.865863 


0.615969 


0.567738 


0.95 


0.137631 


0.627187 


0.0623058 


0.866686 


0.616287 


0.551873 


0.90 


0.137962 


0.614355 


0.0620944 


0.869331 


0.617310 


0.535058 


0.85 


0.138559 


0.603917 


0.0617145 


0.874109 


0.619154 


0.517216 


0.80 


0.139466 


0.596127 


0.0611395 


0.881390 


0.621957 


0.498264 


0.75 


0.140736 


0.591305 


0.0603394 


0.891630 


0.625885 


0.478110 


0.70 


0.142431 


0.589855 


0.0592799 


0.905387 


0.631136 


0.456652 


0.65 


0.144629 


0.592298 


0.0579219 


0.923359 


0.637951 


0.433781 


0.60 


0.147422 


0.599309 


0.0562217 


0.946423 


0.646627 


0.409377 


0.55 


0.150922 


0.611785 


0.0541300 


0.975705 


0.657530 


0.383309 


0.50 


0.155270 


0.630931 


0.0515928 


1.01267 


0.671126 


0.355436 


0.45 


0.160640 


0.658405 


0.0485518 


1.05927 


0.688014 


0.325609 


0.40 


0.167257 


0.696539 


0.0449463 


1.11817 


0.708986 


0.293667 


0.35 


0.175409 


0.748717 


0.0407171 


1.19312 


0.735119 


0.259447 


0.30 


0.185482 


0.820038 


0.0358145 


1.28960 


0.767941 


0.222783 


0.25 


0.198006 


0.918565 


0.0302122 


1.41609 


0.809731 


0.183524 


0.20 


0.213754 


0.105795 


0.0239349 


1.58664 


0.864147 


0.141550 


0.15 


0.233939 


0.126381 


0.0171106 


1.82719 


0.937692 


0.968152 


0.10 


0.260769 


0.159292 


0.0100787 


2.19478 


1.04405 


0.494435 


0.05 


0.299688 


0.222228 


0.00364434 


2.86687 


1.22315 


0.000000 


0.00 


0.500000 


+oo 


0.00000000 


+oo 





Table 2. Values of axis ratios, 631 and £21, and normalized rotation parameters, £ ro t, h,N, vn, 
Xn, xn, related to a sequence of virial equilibrium configurations with anisotropy parameters, 
£11 = C22 = C33 = 1/3. The polar axis ratio of the initial nonrotating configuration is £31 = 1. 
Values at the bifurcation point are listed on the boldface line. 
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Table 3. Values of axis ratios, £31 and £21, and normalized rotation parameters, £ ro t, hN, vn, 
A ^r v»7. related to a sennenee of virial ermilihrinm configurations with anisotronv narame- 
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Table 4. Values of effective anisotropy parameters, £ 33 , polar axis ratios of non rotating con- 
figurations, (e 3 i )i, and polar axis ratios, (e 3 i)b, and normalized rotation parameters, {£ ro t)b, 
(fiN)b, (vN)b, (A]v)fc, (XN)b, of configurations at the bifurcation point, related to sequences 
of virial equilibrium configurations with constant effective anisotropy parameters. Values at 
the isotropic (Cn = (22 = (33) are listed on the boldface line. 
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Table 5. Values of the truncated, scaling radius, S, and some profile parameters, expressed 
in Sect.|2] and particularized to both NFW and MOA density profiles, after comparison with 
the results of high-resolution simulations (Fukushige & Makino 2001). See text for more de- 
tails. The last three parameters, related to the rotational velocity field, have been calculated 
in the special case of constant rotational velocity on the equatorial plane (axisymmetric con- 
figurations) related to either rigid (r] anm — 1r\ rot — 1) or differential (to ensure constant 
rotational velocity everywhere: r/ anm = 3n/8, 7] ro t = 3/4) rotation of the generic, isopycnic 
surface. In the special case of rigidly rotating, triaxial configurations, v anm = v rot = Vi nr , 
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Table Bl. Values of shape factors, A p and (B se i) pp , related to limiting configurations defined 
by the values of the axis ratios, £21 and £31. 



